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Abstract 

We give an exact expression for the distribution of the position X(t) of a single 
second class particle in the asymmetric simple exclusion process (ASEP) where initially 
the second class particle is located at the origin and the first class particles occupy the 
sites Z+ = {1,2,...}. 



1 Introduction 



The asymmetric simple exclusion process (ASEP) |3j is one of the simplest models of 
nonequilibrium statistical mechanics and has been called the "default stochastic model for 
transport phenomena" |8j. A useful concept in exclusion processes is that of a second class 



Imagine that the particles in the system are each called either first class or second 
class. The evolution is the same as before, except that if a second class particle 
attempts to go to a site occupied by a first class particle, it is not allowed to 
do so, while if a first class particle attempts to move to a site occupied by a 
second class particle, the two particles exchange positions. In other words, a first 
class particle has priority over a second class particle. This rule has no effect on 
whether or not a give site is occupied at a given time. The advantage, though, is 
that viewed by itself, the collection of first class particles is Markovian, and has 
the same law as the exclusion process. The collection of second class particles is 
clearly not Markovian. However, the collection of first and second class particles 
is Markovian, and again evolves like an exclusion process. 

Here we consider ASEP on the integer lattice Z with jumps one step to the right with rate 
p and jumps one step to the left with rate q = 1 — p. We assume a leftwards drift, i.e. 
q > p. We further assume that the system has one second class particle initially located at 
the origin and first class particles initially located at sites in 



With the above initial condition, we denote by X{t) the position of the second class particle 
at time t. The purpose of this note is to give an exact expression for the probability that 
the second class particle is at position x at time t, i.e. Py (X(t) = x). (The subscript Y 
denotes the sites of the initial configuration of the first class particles.) Our main result is 
for Y = Z + and is given below in ([9]) and in a slightly different form in (fill . 

J The following quote is taken from Liggett [3]. 




Y = {0 < y x < y 2 < ■ ■ ■ } C Z 
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2 A Basic Lemma 



The single second class particle located at X(t) can be viewed as the (single) discrep- 
ancy under basic coupling between two asymmetric simple exclusion processes rjt and ^ 
where ( t (X(t)) = 1 and rj t (X(t)) = and initially {x : CoO) — 1} — Y' — {0} U Y and 
{x : r] (x) = 1} = Y 

The following identity we first learned from H. Spohn [5] but presumably it has a long 
history: 

Py (X(t) = X)= Py, (Ct(x) = 1) " Py {Vt{x) = 1) . (1) 

For the convenience of the reader, we give a short proof of (pQ) . Let £t and r\t be as above 
evolving together under the basic coupling [2j [3]. Recall that the coupled processes satisfy 
Vt < Ct for all t > since they satisfy this inequality at t = [21 S]@ Define 



J c (x,t) 

X(x,t) 
By counting 
Since 



/]r)t(z) = number of particles in configuration^ with positions < x, 

z<x 

2_,Ct( z ) = number of particles in configuration Ct with positions < x, 



z<x 



1 if X(t) < x, 
if X(i) > x. 

J c (x,t) = J v (x,t)+l(x,t). 

Ey, (J C (x,t)) = ^Ey, (CtO)) = J2 ¥ Y'((t(z) = I) 

z<x z<x 

Ey {J v {x,t)) = J^Ey {r] t {z)) = J^Mvti*) = 1) 

z<x z<x 
Ey (l(x,t)) = Py (X(t) < X) = ^Py(X(t) = z), 



(2) 



z<x 



the expectation of p]) gives 



z<x 



z<x 



z<x 



2 Given two configurations r], ( 6 {0, 1} Z we say r] < ( if r)(x) < ((x) for all i£l 
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from which (pp) follows. 



3 Probability for a site to be occupied in ASEP 



For ASEP with particles initially at Y we denote by x m {t) the position of the mth left- 
most particle at time t (so x m (0) = y m ). In Theorem 5.2 of [6J the authors gave an exact 
expression for Py(x m (i) — x). To state this result we first recall the definition of the r- 
binomial coefficients. For < r := p/q < 1 define for each n G Z + 



n 



1 - T 



In 



]! = [n][n-l] •••[!], [0]! := 1, 



n 



[k)\[n - k)\ 



< k < n, 



and if k > n we set [™] = 0. Equation (5.12) of [6] can be written in the following wa 



Py {X m {t) = X) 



\Y\ 

fc=l scr 

|S|=fe 



(3) 



where, if S := {si, . . . , s^} then 

? fc(fc-l)/2/_-QTO+l r m(m-l)/ 



a(£,Y) 

/(x,fc,0 
e(0 



k-1 

k — in 

#{(s,y) : s e 3,y <E Y, and y < s} 

sum of the positions of the elements of S 1 in Y, 

6-6 ( ii )ir ; 



n 

l<i<j<k 
P 



p + - & 
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3 We make some changes in the notation in (5.12) of [6j. The (p, g)-binomial coefficient [£] of |6] equals 
qk(n-k) times the r-binomial coefficient [£] defined above. The second change is a little more subtle. The 
sum in (5.12) is over all finite subsets S C {1,2,..., \Y\} with \S\ > m. If S — {si, . . . , Sk} the subset 

:= {j/si, ■ • ■ iUs k } and the factor HieS^ - ^ appears in the integrand of (5.12). Thus we can equivalently 
sum over all finite subsets S C Y where now the factor rii<i<fe 4 ^ appears in the integrand. The factor 
a {S) = SieS* °f (5-12) becomes a(S,Y) given above. 

4 All contour integrals are to be given a factor of l/2ni. 
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and Cr is a circle of radius R centered at the origin with R 3> 1 so that all (finite) singularities 
of the integrand are enclosed by Cr. Observe that c m ^ = when m> k. 



Since 



\y\ 



F Y ( Vt (x) = !) = £> (x m (t) = x) 



(4) 



m=l 



we sum the right side of j3]) over all m < k. To carry out this sum recall the r-binomial 
theorem 

n r - 

n 

J. 



E 

j=0 



(_l)jW(j-l)/2 = (1 _ z )(l — ZT )■••(! — ZT 11 - 1 ). 



Using this a simple calculation shows 

k 

y^(_i) m + i r m (' m - i )/ 2 r 



fcm 



m=l 



k — m 



fc-1 



'_l^k+l r -k(k+l)/2 JJ^ _ r j 



Thus 



li'l 



fc-i 



p y (^(x) = 1) = J2(-l) k+1 q k(k - 1)/2 r- k{k+1)/2 r>) 

k=l j=l 

scy JCr JCr i=l 



(5) 



Remark: The above formula holds for either \Y\ finite or infinite. For \Y\ = N, the integral 
of order N in §E§ is gotten from the summand S = Y. Since a(Y, Y) = N(N + l)/2, we get 
for the coefficient of this integral 



JV-l 



N-l 



(6) 



3=1 



4 Probability for a site to be occupied by a second class 
particle 

As above, suppose that our initial configuration consists of a second class particle at site 
and first class particles at sites in Y. As above, set Y' = Y U {0}. The process ( t has 
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initially its particles at sites in Y' . We apply formula (jHj) to the initial configurations Y' and 
Y and by ([I]) we subtract to obtain Py(X(t) = x). If \Y' \ = N there is one iV-dimensional 
integral that comes from the expansion of Py/ = 1) when S = Y' . The coefficient of 
the integral of highest order equals (J@J). 

We now consider the special case of step initial condition] that is, Y = Z + , and use 
Corollary (5.13) of [6] to obtain a more compact expression for Fz+(x m (t) = x). To find 
¥i+{7] t {x) = 1) we again apply (J4|) but use (5.13) of [6j. As above we interchange the sums 
over k and m, use the r-binomial theorem ([4], pg. 26), to conclude 

W> z+ ( Vt (x) = 1) = - TT IB 1 - T ') / • • • / d& • • • dfc (7) 

*>i fc ' i=l ^ 

where 

We can get the corresponding formula for Y 7 = Z + U {0} by observing that there is a one-one 
correspondence between subsets 5" C Y' and subsets S C Y given by S = S' + 1. Then 
a(S', Y') = a(S, Y) and, with obvious notation, riC Si — Yld ' II C - It follows that for 
the difference FyiCti 30 ) = 1) — Py( 7 7t( a; ) — 1) we multiply the integrand Jk(x,£) in by 

lie i. 

Thus 

^>2 fa 1 n n 

P z+ (X(t)=x) = ^^-n( 1 - r ') / •••/ J*(*.£)dei"-d& ( 8 ) 

where 

. p. — p. „ „ t 1 " 1 

y p + q&tj -6 \ x A Mi - 6) (?& - p) 

From this it follows that the distribution function is (on Cr, <C 1) 

P z+ (X(t) < x) = ]T q — - T j ) / • ■ ■ / J k (x, d& ■ ■ ■ dfc (9) 

fc>l K - j=l Jc * Jc * 

where 
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Since 

1 1 

p + q££ - € f (f - i) 

the TASEP limit of J k (x,£) is 

jr Ep (*, o := lim j k ( X , o = n& - (ii & - 1) n 

where now e(£) = £ — 1; and hence, 

limP z+ (X(t) < *) = £ n / j J k ASEF (^0 d£i ■ ■ -d&. (10) 



fc>l 



Expression ((9|) for the distribution function can be simplified somewhat. Define the kernel 

P + - f 

and the associated operator 1^ on L 2 {Cr) by 
Then using the identity [7] 

det f i A r ) = (-i) fe (M) fc(fc - i)/2 n A"/' . n n — T 

we have 

Pz + (X(t) < x) = r- h ^ 1)/2 ;Q(1 - r j ) x 

fc>l j=l 
(-1 



- [ ■ ■ ■ f [det {K x+X ^i, Cj)) 1<iij<k - det Q) l<iyj<k 
Jc R Jc R l 
fc-i 



= ]T r-^*" 1 )/ 2 " / t^t [det(J - A^ +lji ) - det (J - Ai^, t )] rfA (11) 

k >l j=l 

where det(7 — \K X}t ) is the Fredholm determinant and the last line follows from the Fredholm 
expansion. 
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Remarks: 



1. One cannot interchange the sum and the integration in (TTT1) as was possible in an 
analogous calculation in [7]. This is the case even though ([HD converges absolutely for 
all < r < 1 (recall one may take R 1). Thus we do not have a representation 
of F z +(X(t) < x) as a single integral whose integrand involves the above Fredholm 
determinants as was the case in [7]. 

2. ASEP with first and second class particles is integrable in the sense that the Yang- 
Baxter equations are satisfied p]. Using this integrable structure, it is possible to 
compute directly, i.e. without using the basic lemma (pQ), Fz+(X(t) = x) using methods 
similar to that of [6]. We have carried this out to the extent that © was computed 
by this approach. However, this route is much more involved than the one presented 
here. 
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